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Abstract 

This work gives a detailed investigation of matrix product state (MPS) representations for pure 
multipartite quantum states. We determine the freedom in representations with and without transla- 
tion symmetry, derive respective canonical forms and provide efficient methods for obtaining them. 
Results on frustration free Hamiltonians and the generation of MPS are extended, and the use of 
the MPS-representation for classical simulations of quantum systems is discussed. 
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1 Introduction and Overview 

The notorious complexity of quantum many-body systems stems to a large extent from the expo- 
nential growth of the underlying Hilbert space which allows for highly entangled quantum states. 
Whereas this is a blessing for quantum information theory — it facilitates exponential speed-ups in 
quantum simulation and quantum computing — it is often more a curse for condensed matter theory 
where the complexity of such systems make them hardly tractable by classical means. Fortunately, 
physical interactions are local such that states arising for instance as ground states from such inter- 
actions are not uniformly distributed in Hilbert space. Hence, it is desirable to have a representation 
of quantum many-body states whose correlations are generated in a 'local' manner. Despite the fact 
that it is hard to make this picture rigorous, there is indeed a representation which comes close to 
this idea — the matrix product state (MPS) representation. In fact, this representation lies at the 
heart of the power of the density matrix renormalization group (DMRG) method and it is the basis 
for a large number of recent developments in quantum information as well as in condensed matter 
theory. 

This work gives a detailed investigation of the MPS representation with a particular focus on the 
freedom in the representation and on canonical forms. The core of our work is a generalization of 
the results on finitely correlated states in [T | to finite systems with and without translational invari- 
ance. We will mainly discuss exact MPS representations throughout and just briefly review results 
on approximations in Sec(6] In order to provide a more complete picture of the representation and 
its use we will also briefly review and extend various recent results based on MPS, their parent 
Hamiltonians and their generation. The following gives an overview of the article and sketches the 
obtained results: 

• Sec[2]will introduce the basic notions, provide some examples and give an overview over 
the relations between MPS and the valence bond picture on the one hand and frustration free 
Hamiltonians and finitely correlated states on the other. 

• In Sec[3]we will determine the freedom in the MPS representation, derive canonical forms 
and provide efficient ways for obtaining them. Cases with and without translational invari- 
ance are distinguished. In the former cases we show that there is always a translational 
invariant representation and derive a canonical decompositions of states into superpositions 
of 'ergodic' and periodic states (as in ID"). 
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• Sec[4] investigates a standard scheme which constructs for any MPS a local Hamiltonian, 
which has the MPS as exact ground state. We prove uniqueness of the ground state (for 
the generic case) without referring to the thermodynamic limit, discuss degeneracies (spon- 
taneous symmetry breaking) based on the canonical decomposition and review results on 
uniform bounds to the energy gap. 

• In Sec[5] we will review the connections between MPS and sequential generation of mul- 
tipartite entangled states. In particular we will show that MPS of sufficiently small bond 
dimension are feasible to generate in a lab. 

• InSeclU we will review the results that show how MPS efficiently approximate many impor- 
tant states in nature; in particular, ground states of ID local Hamiltonians. We will also show 
how the MPS formalism is crucial to understand the need of a large amount of entanglement 
in a quantum computer in order to have a exponential speed-up with respect to a classical 
one. 



2 Definitions and Preliminaries 

2.1 MPS and the valence bond picture 

We will throughout consider pure quantum states \ip) £ C® d characterizing a system of iV sites 
each of which corresponds to a d-dimensional Hilbert space. A very useful and intuitive way of 
thinking about MPS is the following valence bond construction: consider the N parties ('spins') 
aligned on a ring and assign two virtual spins of dimension D to each of them. Assume that every 
pair of neighboring virtual spins which correspond to different sites are initially in an (unnormal- 
ized) maximally entangled state |J) = J^„ =1 \cx, a) often referred to as entangled bond. Then 
apply a map 

d D 

-4=EE Ai, a ,p\i)(a,/3\ (1) 

1 = 1 a, ,3 = 1 

to each of the iV sites. Here and in the following Greek indices correspond to the virtual systems. 
By writing A; for the D x D matrix with elements Ai !Ctl p we get that the coefficients of the final 
state when expressed in terms of a product basis are given by a matrix product tr [Ai 1 Ai 2 ■ ■ ■ Ai N ] . 
In general the dimension of the entangled state \I) and the map A can both be site-dependent and 
we write A^ for the Dk X Dk+i matrix corresponding to site k 6 {1, . . . , JV}. States obtained in 
this way have then the form 

\i>} = £ tr[A^...Arj]\ il ,i 2 ,...M), (2) 

i 1 ,...,i JV = l 

and are called matrix product states 0. As shown in (7) every state can be represented in this way 
if only the bond dimensions Dk are sufficiently large. Hence, Eq.{2} is a representation of states 
rather than the characterization of a specific class. However, typically states are referred to as MPS 
if they have a MPS-representation with small D — maxfc Dk which (in the case of a sequence 
of states) does in particular not grow with JV. Note that ip in Eq.l(2j is in general not normalized 
and that its MPS representation is not unique. Normalization as well as other expectation values of 
product operators can be obtained from 



(VKXjSfelV) = tr 



with 



4 fc ) = 52(i\S\3)4" ] ®A?K (3) 
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Figure 1 : Computing an expectation value of an MPS is equivalent to contract the tensor of the figure, where 
bonds represent indices that are contracted. The matrices associated to each spin are represented by the circles 
(the vertical bond of each matrix is its physical index) and observables are represented by squares. It is trivial to 
see that this contraction can be done efficiently. 

2.2 Finitely correlated states 

The present work is inspired by the papers on finitely correlated states (FCS) which in turn gener- 
alize the findings of Affleck, Kennedy, Lieb and Tasaki (AKLT) (3). In fact, many of the results 
we derive are extensions of the FCS formalism to finite and/or non-translational invariant systems. 
For this reason we will briefly review the work on FCS. A FCS is a translational invariant state 
on an infinite spin chain which is constructed from a completely positive and trace preserving map 
E : B(Ha) — » B(Ha ® Hb) and a corresponding fixed point density operator A = tr_g[E(A)]. 
Here Hb = C d is the Hilbert space corresponding to one site in the chain and Ha = C D is an 
ancillary system. An n-partite reduced density matrix p n of the FCS is then obtained by repeated 
application of E to the ancillary system (initially in A) followed by tracing out the ancilla, i.e., 

Pn = trA [E n (A)] . (4) 

An important instance are purely generated FCS where E(x) = V^xV is given by a partial 
isometry V. The latter can be easily related to the A's in the matrix product representation via 
V = X/j=i 0=1 Ai lCll p\a>}{f3i\. Expressed in terms of the matrices Ai the isometry condition 
and the fixed point relation read 

d d 

Y,AiA\ = l, £)4AAi=A, (5) 

t=i !=i 

which already anticipates the type of canonical forms for MPS discussed below. As shown in 
141 purely generated FCS are weakly dense within the set of all translational invariant states on 
the infinite spin chain. Moreover, a FCS is ergodic, i.e., an extreme point within all translational 
invariant states, iff the map £(x) = J^, AixA\ has a non-degenerate eigenvalue 1 (i.e., 1 and A 
are the only fixed points in Eq.<[5}). Every FCS has a unique decomposition into such ergodic FCS 
which in turn can be decomposed into p p-periodic states each of which corresponds to a root of 
unity exp(2^pm), m = 0, . . . ,p — 1 in the spectrum of £. A FCS is pure iff it is purely generated 
and 1 is the only eigenvalue of £ of modulus 1. In this case the state is exponentially clustering, 
i.e., the connected two-point correlation functions decay exponentially 

{Si ® I®'" 1 ® S i+ i) - (S t }(S t+ i} = 0(H 1_1 ) , (6) 

where V2 < 1) is the second largest eigenvalue of £. 
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2.3 Frustration free Hamiltonians 



Consider a translational invariant Hamiltonian on a ring of N d-dimensional quantum systems 

N 

H =Y, ri ( h )> (7) 

i=l 

where r is the translation operator with periodic boundary conditions, i.e., r ( &> i=1 Xi) = ® H=i fi+i 
where sites JV + 1 and 1 are identified. The interaction is called L-local if h acts non-trivially only 
on L neighboring sites, and it is said to be. frustration free with respect to its ground state 4>o if the 
latter minimizes the energy locally in the sense that {<j>o \H\ <f>o) — iid^,{4>\H\4>) — N inf^ ((f>\h\(f>) . 
As proven in all gapped Hamiltonians can be approximated by frustration free ones if one allows 
for enlarging the interaction range L up to C(log N). 

For every MPS and FCS ip one can easily find frustration free Hamiltonians such that ip is their 
exact ground state. Moreover, these parent Hamiltonains are L-local with L ~ 2 log D/ log d and 
they allow for a detailed analysis of the ground state degeneracy (Sec l4.lt and the energy gap above 
the ground state (Sec l4.2t . Typically, these Hamiltonians are, however, not exactly solvable, i.e., 
information about the excitations might be hard to obtain. 



2.4 Examples 

1. AKLT: The father of all matrix product states is the ground state of the AKLT-Hamiltonian 

H — ^2 SiSi+i + — (SiSi+i) , (8) 



where S is the vector of spin-1 operators (i.e., d=3). Its MPS representation is given by 
{Ai} = {o -2 , \^2a + , — \f2a~ } where the ct's are the Pauli matrices. 

2. Majumdar-Gosh: The Hamiltonian 

H = 2&iai + i + ffiOi+2 (9) 

i 

is such that every ground state is a superposition of two 2-periodic states given by products of 
singlets on neighboring sites. The equal weight superposition of these states is translational 
invariant and has an MPS representation 
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Ai = -1 , A 2 = 1 . (10) 



3. GHZ states of the form \^>) — | + +...+} + . . . — ) have an MPS representation 

A± = 1 ± cr z . Anti-ferromagnetic GHZ states would correspond to A± — a . 

4. Cluster states are unique ground states of the three-body interactions erf <7f +1 of +2 and 
represented by the matrices 

( \ , ( 1 -1 

M = i ) - A2 = {o 

5. W -states can for instance appear as ground states of the ferromagnetic XX model with strong 
transversal magnetic field. A W-state is an equal superposition of all translates of 1 100 . . . 00} . 
For a simple MPS representation choose {A^\ A^} equal to {cr + , 1} for all k < N and 
{a + a x , a x } for k = N. Although the state itself is translational invariant there is no MPS 
representation with D — 2 having this symmetry. 
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3 The canonical form 



The general aim of this section will be to answer the following questions about the MPS represen- 
tation of a given pure state: 

Question 1 Which is the freedom in the representation? 

Question 2 Is there any canonical representation ? 

Question 3 If so, how to get it? 

We will distinguish two cases. The general case, or the case of open boundary conditions (OBC) 
and the case in which one has the additional properties of translational invariance (TI) and periodic 
boundary conditions (PBC). 

3.1 Open boundary conditions 

A MPS is said to be written with open boundary conditions (OBC) if the first and last matrices are 
vectors, that is, if it has the form 

W)= E ^...A^A^-iN), (ID 

il,...,i N 

where A' m ' are D m x D m +\ matrices with D\ — Dn+i = 1- Moreover, if D — max m D m we 
say that the MPS has (bond) dimension D. The following is shown in 171 : 

Theorem 1 (Completeness and canonical form) Any state ip £ <C d ® N has an OBC-MPS repre- 
sentation of the form Eq. Ml\ with bond dimension D < d'- JV / 2 J aru f 

1. Ei4 ml 4 mlt = lD m forall l<m<N. 

2 ^ A [™]t A [m-i] A [m] = h\. m \forall 1 < m < N, 

3. A [0] = A [N] = 1 and each A [m] is a Djn+i x D m +i diagonal matrix which is positive, full 
rank and with trA [m l = 1. 

ThmCQis proven by successive singular value decompositions (SVD), i.e., Schmidt decompositions 
in ip, and the gauge conditions 1.-3. can be imposed by exploiting the simple observation that 
A M A ljn+n = (A [ ^ n] X)(X- 1 A l J n+1] ). If 1.-3. are satisfied for a MPS representation, then we 
say that the MPS with OBC is in the canonical form. From the way it has been obtained one 
immediately sees that: 

• it is unique (up to permutations and degeneracies in the Schmidt Decomposition), 

• A' m ' is the diagonal matrix of the non-zero eigenvalues of the reduced density operator p m = 

tr m+ i,...,jv W(ip\, 

• any state for which max m rank(p m ) < D can be written as a MPS of bond dimension D. 

This answers questions|2]and[3] Question[T]will be answered with the next theorem which shows 
that the entire freedom in any OBC-MPS representation is given by 'local' matrix multiplications. 

Theorem 2 (Freedom in the choice of the matrices) Let us take a OBC-MPS representation 

Then, there exist (in general non-square) matrices Yj, Zj with YjZj = 1 such that, if we define 
A^=B^Z U AW=Y N ^BW 

4 ml = Y m ^Bl m] Z m , for l<m<N (12) 
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the canonical form is given by 

\*)= E d3) 

ix,...,i N 

Proof. We will prove the theorem in three steps. 

STEP 1. First we will find the matrices Ap' verifying relation d 12b but just with the property 

To this end we start from the right by doing SVD: B [ "} = J2p ^/T^ A Jf '^^l with 
jjIn-D^In] unitariesand A [iv-i] diagonaL ThatisBf 1 =Z N - 1 A\ N] ,withZ N - 1 = U^-n^-n 
Clearly £\ A^A™ = 1 and Z N -x has a left inverse. Now we call b\ n ~ 1] = B^'^Zn^ 
and make another SVD: B 1 "^ = E 7 U^y 21 A^* 21 A^" 1] . That is 

where Af "^Af" 1 " = 1 and Ziv_ 2 = U lN ~ 2] A^" 2 ' has left inverse. 

We can go on getting relations J 12b to the last step, where one simply defines = s] 1 ' Zi . 
From the construction one gets Eq.J13t and that £\ Ap' A^ = 1 for every 1 < j < N, The case 
j = 1 comes simply from the normalization of the state: 

i = <v#>= E AW...<]<^...4 1 1 It =E<< t ' 

where in the last equality we have used that A| J ''A^ = 1 for 1 < j < N. 

STEP2. Nowwecan assume thatthe B'sverify J] i Bp 11 = 1. Diagonalizing J2i B t ^ B f ] 
we get a unitary V [1] and a positive diagonal matrix A [1] such that £\ b\ 1]J< B\ 1] = l/ [11 A [l1 V [llt . 
Calling A- 1 ' = Bf 1 V [1] we have both £\ A™ A™ 1 " = 1 and £\ A^A™ = A' 1 ). 

NowwediagonalizeEiBf'V'^AwVwtBf 1 = V^A^V^ and define A® = V^B^V^ 
to have both A| 2 U| 2]t = 1 and £\ Af ]t A [1 Uf ] = A [2] . We keep on with this procedure to 
the very last step where we simply define A\ N ^ = V^ -1 ^ B\ N \ Aj 4j = 1 is trivially 
verified and £\ A[ m A lN ~ 1] A 1 ^ 1 = A^ = 1 comes, as above, from the normalization of the 
state. Moreover, by construction we have the relation J12b and Eq.J13l>. 

STEP 3. At this point we have matrices Yj, Zj with YjZj = 1 such that, if we define A^' by 
d 1 2b . we get Dj x Dj+i matrices verifying the conditions 1, 2 and 3 of TheoremQ]with the possible 
exception that the matrices A^' are not full rank. Now we will show that we can redefine Yj, Zj 
(and hence Dj, to guarantee also this full rank condition. 

We do it by induction. Let us assume that A'-' -1 ' is full rank and the Dj+i x Dj positive 
diagonal matrix 

A bi 

is not. Then, calling 

D J+1 = rank(A W ) , P 3 = (l 5 . +i \0 D . +1 _ 5j+1 ) , 

we are finished if we update Zj as ZjPj , Yj as PjYj (and hence -Dj+i as Dj+x, A^ as A^Pj, 
A [ * +1] as PjAf +1] and A [jl as PjA^Pj), The only non-trivial part is to prove that A^A^** 1 = 
A i ]p j P 3 A i +1 i' For tnat - calling C = 1d ]+1 - Pj Pj,h is enough to show that A^C = 0. Since 

l 5 

In n 



i B , +1 -p;p, 

we have 



= CA [J, C = ^ CA[ 3lt A [j - 1] Afc. 
i 

Since A' J_1 ' is positive and full rank we get A[^C — 0. q 
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3.2 Periodic boundary conditions and translational invariance 

Clearly, if the A's in the MPS in Eq.© are the same, i.e., site-independent (A^ = Ai), then the 
state is translationally invariant (TI) with periodic boundary conditions (PBC). We will in the fol- 
lowing first show that the converse is also true, i.e., that every TI state has a TI MPS representation. 
Then we will derive canonical forms having this symmetry, discuss their properties and show how 
to obtain them. An important point along these lines will be a canonical decomposition of TI states 
into superpositions of TI MPS states which may in turn be written as superpositions of periodic 
states. This decomposition closely follows the ideas of [T] and will later, when constructing parent 
Hamiltonians, give rise to discrete symmetry-breaking. 



3.2.1 Site independent matrices 

Before starting with the questions Q~| [2] and [3] we will see that we can use TI and PBC to assume 
the matrices in the MPS representation to be site independent. That is, if the state is TI, then there 
is also a TI representation as MPS0 

Theorem 3 (Site-independent matrices) Every TI pure state with PBC on a finite chain has a 
MPS representation with site-independent matrices A^ — Ai, i.e., 



— ^ tr(Aij • • • Ai f 



in) 



(14) 



If we start from an OBC MPS representation, to get site-independent matrices one has ( in general) 
to increase the bond dimension from D to ND (note the N -dependence). 



Proof. We start with an OBC representation of the state with site-dependent A 
the matrices (for 0<i<d— 1) 



[m] 



and consider 



/ 



This leads to 



Bi = N~" 



d-i 

E 



AN 




A 



I-'] 



A 



[N-l] 



tr(B h ---Bi 



jv-i 



r=0 
d-1 



E 



■ A 



IN] 

iN+j 



M. 



□ 



3=0 i 1 ,...,i JV =0 

where ij = ij-N if j > N. Due to TI of tp this yields exactly Eq. J 14t . 

To explicitly show the iV-dependence of the above construction we consider the particular case 
of the VK-state |10 ... 0) + |01 ... 0) + • • • + |0 . . . 01). In this case the minimal bond dimension 
is 2 as a MPS with OBC. However, if we want site-independent matrices, it is not difficult to show 
that one needs bigger matrices. In fact, we conjecture that the size of the matrices has to grow with 
N (Appendix lAl. 

From now on we suppose that we are dealing with a MPS of the form in Eq.dl4t with the 
matrices Ai of size D x D. In cases where we want to emphasize the site-independence of the 
matrices, we say the state is TI represented or simply a TI MPS. 



*In a similar way other symmetries can be restored in the representation. For example if the state is reflection symmetric 
then we can rind a representation with Ai = and if it is real in some basis then we can choose one with real A j . Both 
representations are easily obtained by doubling the bond dimension D. For the encoding of other symmetries in the A's we 
refer to fT1l6ll24l. 
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3.2.2 MPS and CP maps 



There is a close relation (and we will repeatedly use it) between a TI MPS and the completely 
positive map £ acting on the space of D x D matrices given by 

£{X) = Y J A i XA\. (15) 

i 

One can always assume without loss of generality that the cp map £ has spectral radius equal to 1 
which implies by (8] Theorem 2.5] that £ has a positive fixed point. As in the FCS case stated in 
Eq.(|6j the second largest eigenvalue of £ determines the correlation length of the state and as we 
will see below the eigenvalues of magnitude one are closely related to the terms in the canonical 
decomposition of the state. Note that £ and E% = E; Ai (g> Ai have the same spectrum as they are 
related via 

(/3i|£(|Qi)(a 2 |)|/3 2 ) = </3i,/3 2 |£i|ai,a2> . (16) 

Since the Kraus operators of the cp map £ are uniquely determined up to unitaries, it implies 
that £ uniquely determines the MPS up to local unitaries in the physical system. This is used in 
l9l to find the fixed points of a renormalization group procedure on quantum states. There it is 
made explicit in the case of qubits, where a complete classification of the cp-maps is known. To be 
able to characterize the fixed points in the general case one has to find the reverse relation between 
MPS and cp maps. That is, given a MPS, which are the possible £ that can arise from different 
matrices in the MPS representation? It is clear that a complete solution to question Q] will give 
us the answer. However, though we will below provide the answer in the generic case, this is far 
from being completely general. As a simple example of how different the cp-maps £ can be for the 
same MPS, let us take an arbitrary cp-map £ (X ) = E AiXA\ and consider the associate MPS 
for the case of 2 particles: \tp) — 4 tr{Ai 1 Ai 2 )\i\i2) ■ Now translational invariance means 
permutational invariance and hence it is not difficult to show that there exist diagonal matrices Di 
such that \tp) = J3 i i2 tr(Di 1 Di 2 )\iii2) . This defines a new cp-map £(X) = E 4 DiXD\ with 
diagonal Kraus operators, for which e.g many of the additivity conjectures are true 1101 . 



3.2.3 The canonical representation 

In this section we will show that one can always decompose the matrices of a TI-MPS to a canonical 
form. Subsequently we will discuss a generic condition based on which the next section will answer 
question|2]concerning the uniqueness of the canonical form. 

Theorem 4 (TI canonical form) Given a TI state on a finite ring, we can always decompose the 
matrices Ai of any of its TI MPS representations as 

\ 
X 2 A 2 i j , 

where 1 > A j > for every j and the matrices A\ in each block verify the conditions: 

1. E^Mf = i- 

2. Ej A^ A{ = A j ,for some diagonal positive and full-rank matrices A J . 

3. 1 is the only fixed point of the operator £j (X) = Ej A\XA^. 

If we start with a TI MPS representation with bond dimension D, the bond dimension of the 
above canonical form is < D. 
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Proof. We assume w.l.o.g. that the spectral radius of £ is 1 (this is where the Xj appear) and 
we denote by X a positive fixed point of £ . If X is invertible, then calling Bi — X^zAiX? we 
have ^iBiBj = 1 and hence condition 1. 

If X is not invertible and we write X = E a X a \a){a\, and we call Pr the projection onto the 
subspace R spanned by the \a)'s, then we have that AiP R = P R AiP R for every i. To see this, it 
is enough to show that Ai|a) 6 R for every i, |a). If this does not happen for some j, \/3), then 
j: a X a \a)(a\ - XpAAp)(p\A) t 0. But, since E« X a \a){a\ = E, E a X a Ai\a){a\At we 
have obtained that 

X a Ai\a)(a\Al t 0, 

which is the desired contradiction. 

If we call _R X the orthogonal subspace of R, we can decompose our state as 

\i>) = tr R (A il ---Ai N )\i 1 ---i N } + 

+ Yl ten±(Ai 1 ---A iN )\i 1 ---i N ). 

On the one hand tr R (Ai 1 ■ ■ ■ Ai N ) is given by 

tr(P R A n ■ ■ ■ A tN P R ) = tr(P R A n P R ■ ■ ■ P R A lN P R ) 

which corresponds to a MPS with matrices Bi — P R AiP R of size dim(_R) x dim(i?). On the 
other hand 

tr R ± (Ai 1 ■ ■ ■ Ai N ) — tr(P R ±Ai 1 ---Ai N P R i_) 

= ti(P a x Ai 1 P R ± ■ ■ ■ P R ± A lN P R± ) 

since AiP R ± = P R Aj N P R ± + P R ± Ai N P R ± and the P R in the first summand goes through all 
the matrices Ai j to finally cancel with P R ± . Then we have also matrices d = P R ± Ai N P R ± of 
size dim(_R x ) x dim(7? ) such that we can write out original state with the following D x D 
matrices 

Bi 
d 

For each one of these blocks we reason similarly and we end up with block-shaped matrices 
with the property that each block satisfies 1 in the Theorem. Let us now assume that for one of the 
blocks, the map X h- ► E- BiXBj has a fixed point X / 1. We can suppose X self-adjoint and 
then diagonalize it X — E Q with Ai > • • • > X n . Obviously, 1 — -j^-X is a positive 

fixed point that is not full rank, and this allows us (reasoning as above) to decompose further the 
block Bi in subblocks until finally every block satisfies both properties 1 and 3 in the Theorem. By 
the same arguments we can ensure that the only fixed point of the dual map X i— > E; BjXBi of 
each block is also positive and full rank, and so, by choosing an adequate unitary U and changing 
B l to UBiU\ we can diagonalize this fixed point to make it a diagonal positive full-rank matrix 
A, which finishes the proof of the Theorem. [— j 



Note that Thm[4]gives rise to a decomposition of the state into a superposition of TI MPS each 
of which has only one block in its canonical form and a respective cp-map £ j with a non-degenerate 
eigenvalue 1 (due to the uniqueness of the fixed point). The following argument shows that in cases 
where £j has other eigenvalues of magnitude one further decomposition into a superposition of 
periodic states is possible. 

Examples of states with such periodic decompositions (for p = 2) are the anti-ferromagnetic 
GHZ state and the Majumdar-Gosh state. 
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Theorem 5 (Periodic decomposition) Consider any TI state ip £ C d ® N which has only one block 
in its canonical TI MPS representation (Thm^ with respective D X D matrices {Ai}. If£(X) — 
y"] ,. AiXA] has p eigenvalues of modulus one, then ifp is a factor of N the state can be written as 
a superposition ofp p-periodic states each of which has a MPS representation with bond dimension 
D. Ifp is no factor, then ip — 0. 

Proof. The theorem is a consequence of the spectral properties of the cp map £, which 
were proven in (T). There it is shown that if the identity is the only fixed point of £ , then there 
exists a p £ N such that {cu k }k=i... P with ui — exp — are all eigenvalues of £ with modulus 1. 
Moreover, there is a unitary U = X/fc=i u ' Pk, where {Pk} is a set of orthogonal projectors with 
J2h Pk = 1 such that £{XP k ) = £{X)P k -i for all D x D matrixes X (and cyclic index k). It is 
straightforward to show that the latter implies that 

Vj.fc: A j P k = P k - 1 A j . (17) 

Exploiting this together with the decomposition of the trace tr[. . .] = 2~2k tr[Pk . . . Pk] leads to a 
decomposition of the state \tp) = 53? =1 l^fe) where each of the states \ipk) in the superposition 
has a MPS representation with site-dependent matrices Ai . = P k +j-iAi Pk+j- Hence, each \ijjk) 
is p-periodic and, since PkPi = Sk,iPk, non-zero only if p is a factor of N. r- 1 



3.2.4 Generic cases 

Before proceeding we have to introduce two generic conditions on which many of the following 
results are based on. The first condition is related to injectivity of the map 

d 

T L :X^ Y, te[XA il ---Ai L ]\i 1 ...i L }. (18) 

Note that Fl is injective iff the set of matrices {Ai 1 ■ ■ ■ Ai L : 1 < ii, . . . , < d} spans the entire 
space of D x D matrices. Moreover, if J]\ AiA\ = 1 then evidently injectivity of Tl implies 
injectivity of T L i for all L' > L. To see the relation to 'generic' cases consider d randomly chosen 
matrices Ai. The dimension of the span of their products Ai ± ■ ■ ■ Ai L is expected to grow as d L 
up to the point where it reaches D 2 . That is, for generic cases we expect to have injectivity for 
L > 2 1 ? ■ This intuition can easily be verified numerically and rigorously proven at least for 
d — D — 2. However, in order not to rely on the vague notion of 'generic' cases we introduce the 
following: 

Condition CI: There is a finite number Lq £ N such that Tl is injective. 

We continue by deriving some of the implications of condition CI on the TI canonical form: 

Proposition 1 Consider a TI state represented in canonical MPS form (Thm^. If condition 
CI is satisfied for Lo < TV, then 

1. we have only one block in the canonical representation. 

2. if we divide the chain in two blocks of consecutive spins [1 R], [R+ 1 . . . N], both of them 

with at least Lq spins, then the rank of the reduced density operator p\i..m is exactly D . 

Proof. The first assert is evident, since any X which has only entries in the off-diagonal blocks 
would lead to Tl (X) — 0. To see the other implication we take our D-MPS 

\ip)= tr(A h ---A iN )\i 1 ---i N ), 

ix,...,i N 

and introduce a resolution of the identity 

D 

12 (a\A il ---A iR \/3){l3\A iR+1 ---A iN \a)\i 1 ---i N ) = 

ex., (3 = 1 i\ 
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= ^2 \^c,p)\^a,0), where 

a, 0=1 

\® a ,fi) = {a\A n --- A tR \f3)\i 1 ---i R ), 

i 1} ...,i R 

|*«,/3> = J2 {P\A lR ---A %N \a)\i R+ il---i N ). 

i lt ...,i R 

It is then sufficient to prove that both {|$ aj /3}} and {l^a.fi}} are sets of linearly independent vec- 
tors. But this is a consequence of C 1 : Let us take complex numbers c a ./3 such that ~^2 a pC a ^\ & a ,0) '- 
(the same reasoning for the j'I'a.a)). This is exactly 

r R 5> a!i3 |/3><a| =0. 

a,0 

By CI we have that ~}2 a „ c a ,p\{3)(a\ — and hence c a ,p = for every a, (3. q 

Now we will introduce a second condition for which we assume w.l.o.g. the spectral radius of 
£ to be one: 

Condition C2: The map £ has only one eigenvalue of magnitude one. 

Again this is satisfied for 'generic' cases as the set of cp maps with eigenvalues which are degen- 
erated in magnitude is certainly of measure zero. It is shown in (TJ that this condition is essentially 
equivalent (Appendix^ to condition CI. In particular C2 also implies that there is just one block in 
the TI canonical representation (Thm[4j of the MPS \tp) = i N ^{Ai 1 ' ' ' Ai N )\ii ■ ■ ■ ijv). 

Moreover, condition C2 implies that, for sufficiently large N, we can approximate (which 
corresponds to £ N via Eg. J 1 6ft') by |_R)(L|; where \R) corresponds to the fixed point of £ (that is, 
the identity), and {L\ to the fixed point A of the dual map. 

Introducing a resolution of the identity as above, we have that = ^ a ^ \^a,0) \^0,a), with 
l*«,/s) = E n i jy (a\Ai, ■ ■ ■ At N jii • • ■ i n ). But now 

(*«, / s|*a',/j') = {a'\£%(\/3'}{f3\)\a) = W^W. 

up to corrections of the order \vi \ N ' 2 (where V2 is the second largest eigenvalue of £). This implies 
that with increasing N 



becomes the Schmidt decomposition associated to half of the chain. Hence we have proved the 
following. 

Theorem 6 (Interpretation of A) Consider a TI MPS state. In the generic case (condition C2), 
the eigenvalues of its reduced density operator with respect to half of the chain converge with 
increasing N to the diagonal matrix A ® A with A from the TI canonical form (Thm\4$. 

3.2.5 Uniqueness 

We will prove in this section that the TI canonical form in Thm[4]is unique in the generic case. 
Theorem 7 (Uniqueness of the canonical form) Let 

d-l 

|V>>= *r(B n ---B lN )\i 1 ---i N ) 
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be a TI canonical D-MPS such that (i) condition CI holds, (ii) the OBC canonical representation 
of\ip) is unique, and (Hi) N > 2Lo + D 4 (a condition polynomial in D). Then, if\ip) admits 
another TI canonical D-MPS representation 

d-l 

|V>)= E tK<V--Gw)l*i ••■*">> 

there exists a unitary matrix U such that Bi — UdU for every i ( which implies in the case where 
A is non-degenerate that Bi = d up to permutations and phases). 

To prove it we need a pair of lemmas. 

Proposition 1 Let T, S be linear maps defined on the same vector spaces and suppose that there 
exist vectors Y±, . . . ,Y„ such that 

• T(Y k ) = S(Y k+1 ) for every 1 < k < n - 1, 

• Y\, . . . ,Y n -i are linearly independent, 

• Y n = YJkZl ^kY k . 

Consider a solution x 7^ of the equation AiS n_1 + ■ ■ ■ + \ n ^\X = 1, and define 

I12 = X\x 2 + X 2 x 

fin-l = XiX n ~ 1 + ■ ■ ■ + X„-iX = 1. 

Then, ifY = J2^Zl VkY k , we have that Y # and T(Y) = ±S(Y). 
Proof. Clearly 

n—l n—l 

T (Y) = E ^ T ( Yk ) = E ^S(Y k+1 ) = 

fc=l k=l 
/n—2 n—l 
= S I E VkY k+ l + E X kY k 
\k=l k=l 

= S {XiYi + (A 2 + fii)Y 2 + ■■■ + (A n -i + ^-2)^-1) , 

and this last expression is exactly -S(Y) by the definition of the p's. Moreover, since jU n _i = 1 
and Yi, . . . ,Y n -i are linearly independent we have that Y 7^ 0. rj 

The following lemma is a consequence of fTTI Theorem 4.4.14]. 

Proposition 2 If B,C are square matrices of the same size n x n, the space of solutions of the 
matrix equation 

W(C® 1) = {B®t)W 

is S <8> M BJ where S is the space of solutions of the equation XC — BX. 

We can prove now Theorem[7] 

Proof. By Proposition [Tj we know that the matrices in the canonical OBC representation 
of are of dimension D 2 x D 2 for any Lq < j < N — Lq (in particular there are at least D 4 of 
such j's). From the TI MPS representation of \ip) we can obtain an alternative OBC representation 
by noticing that 

d-l 

W= E b^(B l2 ^l)---(B lN ^l)bfJ\i 1 ---i N ), 

i± ,. .. ,ijv— 
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Doing the same with the C's we have also 



d-l 



E 41 ] (Ci 2 ®i) 



i N =0 



Using now Theorem [2] and the fact that between Lq and N — Lo both A}? , (Si (8 1) and 
(Cj ® 1) are D 2 x Z) 2 matrices, we can conclude that there exist invertible D 2 x D 2 matrices 
Wi... W D 4 such that W k (d ® 1) = (B» <g> l)W k +i for every 1 < k < £> 4 - 1. 

Now take n such that Wi, ■ ■ ■ , W n -i are linearly independent but W n = X^=i -^fcWfc- Let us 
define x and . . . , fj, n -i as in Lemma[T]and = X^fc=i t^-kWk- By Lemma[T]we have 
and W(d <g> 1) = (-Bi <S> 1)W for every z. Now, Lemma[2]implies that there exist R 7^ such 
that i?Ci — ^BiR for every i. 

We can use now that A = £\ BJAB, to prove that ClR^ARd = j±?R t AR. Since 
the completely positive map X 1— ► CjXCi is trace preserving (and R^AR 7^ 0) one has that 



Now, from 1 = ]T\ CiC|, we obtain that J] 4 BiRR^Bj = Bi? 1 . Since the Bt's have only 
one box (Proposition!]} we conclude that RR^ = 1 so that R is a unitary. [— ] 

3.2.6 Obtaining the canonical form 

In the previous section we have implicitly used the "freedom" that one has in the choice of the 
matrices in the generic case. In this section we will make this explicit (answering question[T} and 
will use it to show how to obtain efficiently the canonical form (answering question^- 

Let us take a TI state \ip) £ c £i ® JV sucn that the rank of all the reduced density operators is 
bounded by D 2 . Clearly it can be stored using a MPS with OBC in Nd D 2 x D 2 matrices. If we 
are in the generic case and this state has a canonical form verifying condition CI, it would be very 
convenient to have a way of obtaining it, since it allows us to store the state using only d D x D 
matrices ! 

In this section we will show how the techniques developed so far allow us to do it by solving an 
independent of system of 0(D S ) quadratic equations with 0(D 4 ) unknowns. 

We will assume that the problem has a solution, that is, the state has a TI canonical form with 
condition CI. We will also assume that we are in the generic case in the sense that the OBC 
canonical form is unique (no degeneracy in the Schmidt Decomposition). Then, the algorithm to 
find it reads as follows: 

We start with the D 2 x D 2 matrices A { ^ 0+1] Af 0+Di] of the OBC canonical form. 

We solve the following system (S) of quadratic equations in the unknows Yj, Zj+i,Bi, j — 
Lo + 1, . . . , Lq + D 4 , i = 1, . . . , d (Yj, Zj are D 2 x D 2 matrices and Bi D x D matrices). 




= 1. 



B i ®l = Y j AfZ L 
YjZj = 1 



E B 4 = 1 



and we have the following. 
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Theorem 8 (Obtaining the TI canonical form) Consider any TI state \ with unique OBC canon- 
ical form and such that the rank of each reduced density operator (of successive spins) is bounded 
byD 2 . 

1. If there is a TI MPS representation verifying condition CI, then the above systems (S) of 
quadratic equations has a solution. 

2. Any solution of(S) gives us a TI D-MPS representation of\ip), that is related to the canonical 
one by unitaries (Ai = U BiU^ ). 

Proof. We have seen in the previous Section that the canonical representation is a solution for 
(S). Now, if Bi's are the solution of (S) and Ai's are the matrices of the canonical representation, 
we have, reasoning as in the proof of Theorem[7] that there exists an R 7^ and an x 7^ such that 
RBi = ±AiR for every i. Using that £\ BiBj = 1, that A = £\ A\kAi and that 1 is the only 
fixed point of X 1— > ^\ AiXA] we can conclude, as in the proof of Theorem[7] that R is unitary. 
□ 



4 Parent Hamiltonians 

This section pretends to extend the results of the seminal paper |jTj to the case of a finite chain. That 
is, we will study when a certain MPS is the unique ground state of certain gapped local hamiltonian. 
However, since we deal with a finite chain, the arguments given in (T) for the "uniqueness" part are 
no longer valid, and we have to find a different approach. As in the previous section we will start 
with the case of OBC and then move to the case of TI and PBC. In the "gap" part we will simply 
sketch the original proof given in |T). 



4.1 Uniqueness 

4.1.1 Uniqueness of the ground state under condition CI in the case of OBC 

Let us take a MPS with OBC given in the canonical form \ ip) ~ 2~2i 1 i ^1*' ' ' ' ^ijv I* 1 
Let us assume that we can group the spins in blocks of consecutive ones in such a way that, in 
the regrouped MPS \%p) = 2~2i ± j ■ ■ ■ • • • i&i), every set of matrices B^} verifies 

condition CI, that is, generates the corresponding space of matrices. If we call hjj+i the projector 
onto the orthogonal subspace of 

{ J2 ti{XB b] B [3+1] ) : X arbitrary}, 

then 

Theorem 9 (Uniqueness with OBC) \ip) is the unique ground state of the local Hamiltonian H = 
Proof. Any ground state \ cj>) of H verifies that hj t j+i ® l\(j>) — for every j, that is 

il,...,iu 

where j + 1) is the set of indices i\ . . . ij-i, ij+2 ■ ■ Am- 

Mixing J19l > for j and j + 1 and using condition CI for B' J+1 ' gives 

yi r[j] _ R L?+2] yj+i 

Using now that £ Sj B^B^ f = 1 and calling V^... J+2) = E Sj ^J+ij+a)^ 1 f we get 
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I*>= E tT ( Y iu,...j+2) B i- ■ ■ ■ • ■ • 

Using the trivial fact that blocking again preserves condition CI one can easily finish the argu- 
ment by induction. We just notice that in the last step one obtains 

\<f>) = J2 XB^...B^\i 1 ...i M ), 

il,...,i M 

where X is just a number that, by normalization, has to be 1, giving \(j>) — \ip) and hence the result. 
□ 

4.1.2 Uniqueness of the ground state under condition CI with TI and PBC 

To obtain the analogue result in the case of TI and PBC one can apply the same argument. However, 
since we do not have any more vectors in the first and last positions, we need to refine the reasoning 
of the last step. Moreover, using the symmetry we have now, one can decrease a bit the interaction 
length of the Hamiltonian, from 2Lq to Lo + 1. 

Let us be a bit more concrete. Given our ring of TV d-dimensional quantum systems, L < N 
and a subspace S of C d , we denote Hs = Yl =i t'O^s). where hs is the projection onto 
S X Q. If we start with a MPS = J2 il i N M^n ' " ' • • • »jv) with property CI, we 

will consider L > Lo and, as before, the subspace G£ (or simply Ql) formed by the elements 
J3 f i tv{XAi 1 ■ ■ ■ Ai L )\ii ■ ■ ■ il). It is clear that Hg L \ip) = and that Hg L is frustration 
free. Moreover, if TV > 2Lo and L > Lq, then 

Theorem 10 (Uniqueness with TI and PBC) \ip) is the only ground state of Hg L . 

Proof. Reasoning as in the case of OBC one can easily see that any ground state \cf>) of Hg L 
is in Qn, that is, has the form \(j>) = ^ s ijv tr(XAi 1 ■ ■ ■ Ai N )\i\ ■ ■ ■ ijy). Since there is no 
distinguished first position, \<f>) can also be written 

10= ti(A il ---A iLo YA iLo+1 ---A iN )\i 1 ---i N ). 

ix,...,i N 

By condition CI, XAi ± ■ ■ ■ Ai Lg = Ai 1 ■ ■ ■ Ai Lg Y for every ii,... , ih Q . But, also by condition 
CI, Ai 1 ■ ■ ■ Ai Lg generates the whole space of D x D matrices. Hence X = Y which, in addition, 
commutes with Ai x ■ ■ ■ Ai Lg for every i\, . . . , . This means that X commutes with every matrix 
and hence X — XI and \<j>) = \ip). rj 

4.1.3 Non-uniqueness in the case of two or more blocks 

In the absence of condition CI in our MPS, we cannot guarantee uniqueness for the ground state of 
the parent Hamiltonian. There are two different properties that can lead to degeneracy. One is the 
existence of a periodic decomposition (Theorem[5}, that can happen even in the case of one block. 
This is the case of the Majumdar-Gosh model (Section [2}. The other property is the existence of 
more than one block in the canonical form (Thm©. As we will see below, this leads to a stronger 
version of degeneracy that is closely related to the number of blocks. In particular, we are going to 
show that whenever we have more than one block, the MPS is never the unique ground state of a 
frustration free local Hamiltonian (Thm ll lb . In addition, there exists one local Hamiltonian which 
has the MPS as ground state and with ground space degeneracy equal to the number of blocks 

^For all the reasonings it is enough to consider h : <C d ® L > Qd®L positive such that ker h = S. We take h = P s 

simplicity. 
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(Thmll2b. That is, a number of blocks greater than one can correspond to a spontaneously broken 
symmetry (6). 

In all this section we need to assume that we have condition CI in each block. For a detailed 
discussion of the reasonability of this hypothesis see Appendix lAl 

Let us take a TI MPS \<j>) with b(> 2) blocks A}, . . . , A\ in the canonical form, size(A I 1 ) > 
■ ■ ■ > size(j4i), and condition CI in each block. Let us call Lo = maxj{L^ J }, which will be 
logarithmic in N for the generic case. Clearly \<f>) = Ej=i l&Ai }> where 

l^) = Af E tT(Ai 1 ---A{ N )\i 1 -.-i N }. 

i lt ...,i N 

Moreover, w.l.o.g. we can assume that the states \4>Ai) are pairwise different. The following 
lemmas will take care of the technical part of the section. 

Proposition 3 Given any D x D matrices C / and X there exist matrices Ri , Si such that 
X = Y2i RiCSi. 

Proof. By the polar decomposition, it is easy to find matrices E, F, G, H (with G, H invert- 
ible) such that ECF = andGXH = J2?=i |*)(«|- Clearly |*)<*1 = E™i ^l 1 )^!^. 

where Pi is the matrix obtained from the identity by permuting the first and the i-th row. So 
Ri = G~ x P l E and S t = FP.H" 1 . n 

Proposition 4 If L > 3(6 - 1)(L + 1), the sum ©J =1 Ql « direct. 

Proof. We group the spins in 3(6 — 1) blocks of at least Lo + 1 spins each and then use 
induction. First the case 6 = 2. 

We assume on the contrary that there exist X, Y 7^ such that 

E HA^XA^Al^i^is) = E HA^Y^A^hhis). 

If we consider now an arbitrary matrix Z, by CI and Lemma [3] there exist complex numbers 
M | ljP | 2 such that Z = E»E il>i2 A^A\XA\^ Calling W = E.E^ A,pl, A t YA l 
we have that £ l3 tr(Z^ 3 )|i 3 ) = E, 3 tr{W A 2 i3 )\i 3 ) . 

This means that Gl +i C Gl +i- Since size(Aj) > size(Af), this implies that size(A} ) = 
size(Af) and that £?£ n +i = Qtn+i- But now, taking the local Hamiltonian H A i = H A i , 

by Theorem[To] both \§ A i ) and 1^42) should be its only ground state; which is the desired contra- 
diction. 

Now the induction step. Let us start with Ejii ^ — 0, where 

<l 0+ i)3^= E te(M 1 W j Aj 2 ---Alj\i 1 ..-i 3b ). 

ii»...i<3i 

We want to prove that W J = for every j. So let us assume the opposite, take j such that W J 7^ 
and call Wj = I lr 12 i <8> h A b+i ® lu 3i ,i j (to?). We have that E,-=i = 0> an£ ^' by me 
induction hypothesis, each Wj — 0. Now 

•1 «3t 
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with X 3 i3 7^ for some 13 (Theorem 1 10b. Moreover, we can use condition CI and Lemma[3]to get 
complex numbers pl i: pl 2 such that 1 =X)iZ)i I ,i 2 ^Pi^A^W 3 A 3 i2 . 

Hence, = i 3b t r (^i 3 ^4 ' ' ' ^i 36 ) 1*4 ' ' ' *3f>) f° r every 13, which implies by CI the 

contradiction X\ 3 =0. rj 

Finally the results, 

Theorem 11 (Degeneracy of the ground space vl) If N > 3(6 — l)(£o + 1) + L and H — 
y^ . t 1 (/i) is arcv translationally invariant frustration free L-local Hamiltonian on our ring of N 
spins that has \ <j>) as a ground state (that is, H\cj>) = 0), then \ 4>Ai ) * s a ^ so a ground state of H for 
every j. 

In particular H has more than one ground state. 

Proof. One has = (/i<g>l)|0) = 1)1^)- Since (ft® ) £ (C d )® i ®^l L , 

we can use Lemma[4]to get the desired conclusion: (h ® 1) \(f> A j } = for every j. [— j 

Theorem 12 (Degeneracy of the ground space v2) There exists a local Hamiltonian H acting on 
L > 3(6— 1)(Lq + 1) + 1 spins such that its ground space is exactlyker(H) — span{|</> A j)}i<j<b. 



Proof. The hamiltonian will be H s with S = ® j G? , and L > 3(6 - 1)(L + 1) + 1. For 
m> L, 



In fact, if |0) G C d ® (0 J 5^ ) n (©^ ) ® C d , we have simultaneously that 
W=£ £ t *(M m+1 Ci A i--- A L)\ii---irn+i) and 

3 ii,...,i m+ i 

iv>>=£ £ M^ ra+ X< •• -wo ■ 

3 >«tn+l 

Lemma[4]and condition CI allows us to identify for every j, i\ , i m +i 

A3 ni _ p,i a 3 

Calling ^ = E n C?X and using that £ n 4f ; /If / = 1, we get A^C^ = , 

which implies that \<j>) £ (Q. Gm+i- 

Then one can easily follow the lines of the proof of Theorem [To] (assuming N > L + Lq) to 
conclude that ker(_ffg) = span{j</> A j }}i<j<b- rj 



4.2 Energy gap 

If the ground state energy is zero (which can always be achieved by a suitable offset), the energy 
gap 7 above the ground space is the largest constant for which 

H 2 > jH . (20) 

If in addition H — £V T % (h) is frustration free and has interaction length I, by taking any 
p > I and grouping the spins in blocks of p, one can define an associated 2-local interaction in the 
regrouped chain by hi t i+\ = H[ P i+i,...,(i+2) P ] ( = Y^lpP ~' T J (h)). The new hamiltonian H 
verifies 

m — 1 

% rap] < ^[l...m] = E rl (^) — 2^[l,...,mp] • 

i=0 
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Moreover, calling P to the projection onto ker(h), there exists a constant 72 P (that is exactly the 
spectral gap of iSh...,2p]) suc h that h > 72p(l — P)- 

Therefore, to study the existence of an energy gap in a local TI frustration free Hamiltonian, 
it is enough to study the case of a nearest neighbor interaction H — £\ r^Pj.i+i) where P is a 
projector. In this situation, Knabe 1 1 21 gave a sufficient condition to assure the existence of a gap, 
namely that the gap e n of H[i.,., n +i] is bigger than — for some n. 

In the particular case of the parent hamiltonian of an MPS, a much more refined argument was 
provided in [I] to prove the existence of a gap under condition CI. The idea reads as follows. 
Clearly H 2 > H + P.i+iP+i ,i+2 + Pi, 4+2/1+1,1+1. By the proof of Theorem [Tol 

Pi,i+i = l[i...p<] ® (1 — Pg 2p ) ® l[p(i+2)+i...JV]- 

Then a technical argument proves that Pi,i+iPi+i,i+2 + P%,i+aPi+i,i+i > — 0(^f)(-f»,i+i + 
Pi+1,1+2), where V2 is the second largest eigenvalue of £ . This gives H 2 > (1 — Oiv^ ))H, which 
concludes the argument. 

5 Generation of MPS 

The MPS formalism is particularly suited for the description of sequential schemes for the gener- 
ation of multipartite states. Consider for instance a chain of spins in a pure product state. Two 
possible sequential ways of preparing a more general state on the spin chain are either to let an 
ancillary particle (the head of a Turing machine) interact sequentially with all the spins or to make 
them interact themselves in a sequential manner: first spin 1 with 2 then 2 with 3 an so on. 

Clearly, many physical setups for the generation of multipartite states are of such sequential 
nature: time-bin photons leaking out of an atom-cavity system, atoms passing a microwave cavity 
or laser pulses propagating through atomic ensembles. We will see in the following that the MPS 
formalism provides the natural language for describing such schemes. This section reviews and 
extends the results obtained in 1131 . A detailed application of the formalism to particular physical 
systems can be found in 1 1 3 1 1 141 . 

5.1 Sequential generation with ancilla 

Consider a spin chain which is initially in a product state ]0} 8iv € TC^ N with Tin — C d and an 
additional ancillary system in the state \ipi) £ Ha — C D . Let J]\ a a-A^j, g\a, *)(/?, 0| be a 
general stochastic operation on TCa ®7~Lb applied to the ancillary system and the fc'th site of the 
chain. This operation could for instance correspond to one branch of a measurement or to a unitary 
interaction in which case £\ A'^A^ = 1. If we let the ancillary system interact sequentially 
with all N sites and afterwards measure the ancilla in the state \<Pf), then the remaining state on 
Tig is (up to normalization) clearly given by the MPS 

d 

|V) = £ (< PF \A^...AfjA^ I }\i N ...i 1 ). (21) 

ix,...,i N =l 

By imposing different constraints on the allowed operations (and thus on the A's) we can distinguish 
the following types of sequential generation schemes for pure multipartite states: 

1. Probabilistic schemes: arbitrary stochastic operations with a £>-dimensional ancilla are al- 
lowed. 

2. Deterministic schemes: the interactions must be unitary and the D-dimensional ancilla must 
decouple in the last step (without measurement). 
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3. Deterministic transition schemes (for d = 2): Here we consider an enlarged ancillary system 
TLa = C D © C D ~ C c <S> C 2 (corresponding e.g. to D 'excited' and D 'ground state' 
levels) and a fixed interaction of the form 

|<P>|1>|0>» >-» lv>|0>|l>a, 
^)|0)|0> B ^ |V>|0>|0>». 

Moreover, we allow for arbitrary unitaries on TCa in every step and require the ancilla to 
decouple in the last step. 

Theorem 13 (Sequential generation with ancilla) The three sets of multipartite states which can 
be generated by the above sequential schemes with ancilla are all equal to the set of states with 
OBC MPS representation with maximal bond dimension D. 

Note that the proof of this statement in [13] (based on subsequent singular value decompositions) 
also provides a recipe for the generation of any given state (with minimal resources). This idea has 
been recently exploited in 1151 to analyze the resources needed for sequential quantum cloning. 

5.2 Sequential generation without ancilla 

Let us now consider sequential generation schemes without ancilla. The initial state is again a 
product 10}®^ £ c d ® N an( j we perform first an operation affecting the sites 1 and 2, then 2 and 3 
up to one between iV — 1 and N. Again we may distinguish between probabilistic and deterministic 
schemes and as before both classes coincide with a certain set of MPS: 

Theorem 14 (Sequential generation without ancilla) The sets of pure states which can be gen- 
erated by a sequential scheme without ancilla either deterministically or probabilistically are both 
equal to the class of states having an OBC MPS representation with bond dimension D < d. 

Proof. Let us denote the map acting on site k and k + 1 by U^ k \ Then for k < N — 1 we 
can straight forward identify the matrices in the MPS representation by A^ „ = (i, /3|[/' fc ' \a, 0) 
where for k = 1 we have a — 0, i.e., are vectors. From the last map with coefficients 
(i, j\U^ N ~^\a, 0} we obtain the A's by a singular value decomposition (in i\a, j), such that A\ N ^ 
is again a set of vectors. Hence, all states generated in this way are OBC MPS with D < d. 

Conversely, we can generate every such MPS deterministically in a sequential manner without 
ancilla. To do this we exploit Thm[T3]and use the site k + 1 as 'ancilla' for the fc'th step (i.e., the 
application of a unitary C/' fe ') followed by a swap between site k + 1 and k + 2. N — 1 of these 
steps are sufficient since the last step in the proof for the deterministic part in Thmll3lis just a swap 
between the ancilla and site N. i—i 



6 Classical simulation of quantum systems 

We saw in Section|2]that if a quantum many-body state has a MPS representation with sufficiently 
small bond dimension D, then we can efficiently store it on a classical computer and calculate 
expectation values. The practical relevance of MPS representations in the context of classical 
simulation of quantum systems stems then from two facts: (i) Many of the states arising in con- 
densed matter theory (of one-dimensional systems) or quantum information theory either have such 
a small-D MPS representation or are well approximated by one; and (ii) One can efficiently obtain 
such approximating MPS. 

The following section will briefly review the most important results obtained along these lines. 
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6.1 Properties of ground states of spin chains 



The main motivation for introducing the class of MPS was to find a class of wavefunctions that 
capture the physics needed to describe the low-energy sector of local quantum Hamiltonians. Once 
such a class is identified and expectation values of all states in the class can be computed efficiently, 
it is possible to use the corresponding states in a variational method. The very successful renormal- 
ization group methods, first developed by Wilson 1 17) and later refined by White 1181 . are precisely 
such variational methods within the class of matrix product states I19II20II2D . 

In the case of 1-D systems (i.e. spin chains) with local interactions, the low-energy states indeed 
exhibit some remarkable properties. First, ground states have by definition extremal local properties 
(as they minimize the energy), and hence their local properties determine their global ones. Let us 
consider any local Hamiltonian of TV spins that exhibits the property that there is a unique ground 
state \ip ex ) and that the gap is A(7V). Let us furthermore consider the case when A(TV) decays 
not faster than an inverse polynomial in TV (this condition is satisfied for all gapped systems and 
for all known critical translationally invariant systems). Then let us assume that there exists a state 
\ipap P r) that reproduces well the local properties of all nearest neighbor reduced density operators: 
\\pap P r — Pei || < 5. Then it follows that the global overlap is bounded by 

Hl^es) - \^appr)f < ■ 

This is remarkable as it shows that it is enough to reproduce the local properties well to guarantee 
that also the global properties are reproduced accurately: for a constant global accuracy e, it is 
enough to reproduce the local properties well to an accuracy S that scales as an inverse polynomial 
in the number of spins. This is very relevant in the context of variational simulation methods: if the 
energy is well reproduced and if the computational effort to get a better accuracy in the energy only 
scales polynomially in the number of spins, then a scalable numerical method can be constructed 
that reproduces all global properties well (here scalable means essentially a polynomial method)0. 

Second, there is very few entanglement present in ground states of spin chains, even in the 
case of a critical system. The relevant quantity here is to study area-laws: if one considers the 
reduced density operator pl of a contiguous block of L spins in the ground state of a spin chain 
with TV 2> L spins, how does the entropy of that block scale with L? This question was first studied 
in the context of black-hole entropy 1221 and has recently attracted a lot of attention [23 24 25 1. 
Ground states of local Hamiltonians of spins seem to have the property that the entropy is not an 
extensive property but that the leading term in the entropy only scales as the boundary of the block 
(hence the name area-law), which means a constant in the case of a 1-D system I25II26I : 

^( Pi )~£fl + lW(fl (22) 

Here S a is the Renyi entropy 

S a (p) = r 1 - log (Trp a ), 
1 — a 

c is the central charg^land £ the correlation length. This has a very interesting physical meaning: it 
shows that most of the entanglement must be concentrated around the boundary, and therefore there 
is much less entanglement than would be present in a random quantum state (where the entropy 
would be extensive and scale like L). The area law J22l l is mildly violated in the case of 1-D critical 
spin systems where £ has to be replaced with L, but even in that case the amount of entanglement is 
still exponentially smaller than the amount present in a random state. This is very encouraging, as 

t(I)f course this does not apply to global quantities, like entropy, where one needs exponential accuracy in order to have 
closeness 

§ We note that Eq.(22} has been proven for critical spin chains (in particular the XX-model with transverse magnetical field 
1261 ) which are related to a conformal field theory. A general result is still lacking. 
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one may exploit the lack of entanglement to simulate these systems classically. Indeed, we already 
proved that MPS obey the same property 1241 . 

The existence of an area law for the scaling of entropy is intimately connected to the fact that 
typical quantum spin systems exhibit a finite correlation length. In fact, it has been recently proven 
1271 that all connected correlation functions between two blocks in a gapped system have to de- 
cay exponentially as a function of the distance of the blocks. Let us therefore consider a 1-D 
gapped quantum spin system with correlation length £ C orr. Due to the finite correlation length, it 
is expected that the reduced density operator pAB obtained when tracing out a block C of length 
Lab > £,corr is equal to 

Pab ~ pa® Pb (23) 

up to exponentially small corrections The original ground state \tPabc) is a purification of this 
mixed state, but it is of course also possible to find another purification of the form | ipAC t } <8> | ipBC-r ) 
(up to exponentially small corrections) with no correlations whatsoever between A and B; here C'i 
and CV together span the original block C. Since different purifications are unitarily equivalent, 
there exists a unitary operation Uc on the block C that completely disentangles the left from the 
right part: 

Ia®Uc ® Ic\i>ABc) =± \lpACt) ® \lpBCr)- 

This implies that there exists a tensor Ai >a ,p with indices 1 < i,a,/3 < D (where D is the 
dimension of the Hilbert space of C) and states \ipp) defined on the Hilbert spaces 

belonging to A, B, C such that 

\lpABc) ^ ^.«./3|V'£>IV'f>IV , f>- 

Applying this argument recursively leads to a matrix product state description of the state and gives 
a strong hint that ground states of gapped Hamiltonians are well represented by MPS. It turns out 
that this is even true for critical systems. 



6.2 MPS as a class of variational wavefunctions 

Here we will review the main results of [29], which give analytical bounds for the approximation of 
a state by a MPS that justify the choice of MPS as a reasonable class of variational wavefunctions. 
We consider an arbitrary state | ?/>} and denote by 

{p [k]l },i = l..N k = d m ' m(k ' N - k) 

the eigenvalues of the reduced density operators 

Pk = Tr k +i,k+2,...,N|V')(V , | I 

sorted in decreasing order. 

Theorem 15 There exists a MPS \i/jd) with bond dimension D such that 

JV-l 

ll!^)-|^)l| 2 <2^6 fe (D) 
fc=l 

where e k (D)=J2l k D+1 ^ lk]l - 



^Strictly speaking, Hastings theorem does not imply the validity of equation I23K as it was shown in [ 28 ] that orthogonal 
states exist whose correlation functions are exponentially close to each other; although it would be very surprising that ground 
states would exhibit that property, this prohibits to turn the present argument into a rigorous one. 
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This shows that for systems for which the (D) decay fast in D, there exist MPS with small D 
which will not only reproduce well the local correlations (such as energy) but also all the nonlocal 
properties (such as correlation length). 

The next result relates the derived bound to the Renyi entropies of the reduced density operators. 
Given a density matrix p, we denote as before e(D) = X^u+i ^ i w ' m ^> me nonincreasingly 
ordered eigenvalues of p. Then we have 

Theorem 16 I/O < a < 1, then log(e(D)) < 1=2 (s a (p) - log 

The two results together allow us to investigate the computational effort needed to represent 
critical systems, arguable the hardest ones to simulate^], as MPS. The key fact here is the area-law 
i22\ , which for critical systems reads 

S a ( PL )^^l + i)log(L) (24) 

for all a > (c the central charge). 

Let us therefore consider the Hamiltonian associated to a critical system, but restricted to 2L 
sites. The entropy of a half chain (we consider the ground state \ip ex ) of the finite system) will 
typically scale as in eq. d24h but with an extra term that scales like 1/L. Suppose we want to get 
that IH^ei) — IV'd)!! 2 < £o/£ with eo independent of L0. If we call Dl the minimal D needed 
to obtain this precision for a chain of length 2L, the previous two results combined yield 

Theorem 17 

D L <cst[- — LT5- — . 

This shows that D only has to scale polynomially in L to keep the accuracy eo / L fixed; in other 
words, there exists an efficient scalable representation for ground states of critical systems (and 
hence also of noncritical systems) in terms of MPS. This is a very strong result, as it shows that 
one can represent ground states of spin chains with only polynomial effort (as opposed to the 
exponential effort if one would do e.g. exact diagonalization). 

The above result was derived from a certain scaling of the block Renyi entropy with a < 1. 
In fact, for a > 1, i.e., in particular for the von Neumann entropy (a = 1), even a saturating 
block entropy does in general not allow to conclude that states are efficiently approximable by 
MPS 1 16|. Conversely, a linear scaling for the von Neumann entropy, as generated by particular 
time evolutions, rules out such an approximability 1161 . 



6.3 Variational algorithms 

Numerical renormalization group methods have since long been known to be able to simulate spin 
systems, but it is only recently that the underlying structure of matrix product states has been 
exploited. Both NRG, developed by K. Wilson (T7) in the '70s, and DMRG, developed by S. 
White [18) in the '90s, can indeed be reformulated as variational methods within the class of MPS 
12 II 1201 . The main question is how to find the MPS that minimizes the energy for a given spin 
chain Hamiltonian: given a Hamiltonian Tt acting on nearest neighbours, we want to find the MPS 
\ip) such that the energy 

is minimized. If is a TI MPS, then this is a highly complex optimization problem. The main 
trick to turn this problem into a tractable one is to break the translational symmetry and having site- 
dependent matrices A\ for the different spins k; indeed, then the functional to be minimized is a 



'For non-critical systems, the renormalization group flow is expected to increase the Renyi entropies in the UV direction. 
The corresponding fixed point corresponds to a critical system whose entropy thus upper bounds that of the non-critical one. 
"We choose the 1/L dependence such as to assure that the absolute error in extensive observables does not grow. 
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multiquadratic function of all variables, and then one can use the standard technique of alternating 
least squares to do the minimization 1211 . This works both in the case of open and closed boundary 
conditions, and in practice the convergence of the method is excellent. The computational effort of 
this optimization scales as D A d 2 in memory and similarly in time (for a single site optimization). 

But what about the theoretical worst case computational complexity of finding this optimal 
MPS? It has been observed that DMRG converges exponentially fast to the ground state with a 
relaxation time proportional to the inverse of the gap A of the system 1301 . For translationally 
invariant critical systems, this gap seems to close only polynomially. As we have proven that 
D only have to scale polynomially too, the computational effort for finding ground states of 1-D 
quantum systems is polynomial (P). This statement is true under the following conditions: 1) the 
a-entropy of blocks in the exact ground state grow at most logarithmically with the size of the block 
for some a < 1; 2) the gap of the system scales at most polynomially with the system size; 3) given 
a gap that obeys condition 2, there exists an efficient DMRG-like algorithm that converges to the 
global minimum. As the variational MPS approach 1211 is essentially an alternating least squares 
method of solving a non-convex problem, there is a priori no guarantee that it will converge to the 
global optimum, although the occurrence of local minima seems to be unlikely 1301 . But still, this 
has not been proven and the worst-case complexity could well be NP-hard, as multiquadratic cost 
functions have been shown to lead to NP-hard problems 1311 .1^1 

6.4 Classical simulation of quantum circuits 

In the standard circuit model of quantum computation a set of unitary one and two-qubit gates is 
applied to a number of qubits, which are initially in a pure product state and measured separately 
in the end. In the cluster state model a multipartite state is prepared in the beginning and the com- 
putation is performed by applying subsequent single-qubit measurements. For both computational 
models the MPS formalism provides a simple way of understanding why a large amount of en- 
tanglement is crucial for obtaining an exponential speed-up with respect to classical computations. 
The fact that quantum computations (of the mentioned type) which contain too little entanglement 
can be simulated classically is a simple consequence of the following observations (see also 1331 ): 

1. Let £ C® d have a MPS representation with maximal bond dimension D, By Eq.([3} 
the expectation values of factorizing observables are determined by a product of N D 2 x D 2 
matrices. Hence, their calculation as well as the sampling of the respective measurement 
outcomes and the storage of the state requires only polynomial resources in N and D. 

2. A gate applied to two neighboring sites increases the respective bond dimension (Schmidt 
rank) at most by a factor d 2 and the MPS representation can be updated with poly(D) re- 
sources. 

3. A gate applied to two qubits, which are / sites apart, can be replaced by 21 — 1 gates acting 
on adjacent qubits. In this way we can replace a circuit in which each qubit line is crossed 
or involved by at most L two-qubit gates by one where each qubit is affected only by at most 
4L nearest neighbor gates. 

A trivial implication of 1. is that measurement based quantum computation (such as cluster state 
computation) requires more than one spatial dimension: 

Theorem 18 (Simulating ID measurement based computations) Consider a sequence of states 
with increasing particle number N for which the MPS representation has maximal bond dimension 
D = 0(poly(N)). Then every measurement based computation on these states can be simulated 
classically with 0(poly(N)) resources. 

ttln fact, it can be shown |32| that a variant of DMRG leads to NP-hard instances in intermediate steps. However, one has 
to note that (i) such a worst case instance might be avoided by starting from a different initial point and (ii) convergence to the 
optimum at the end of the day does not necessarily require finding the optimum in every intermediate step. 
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This is in particular true for a one-dimensional cluster state, for which D — 2 is independent of N. 
Similarly, if the initial states are build by a constant number of nearest neighbor interactions on a 
two-dimensional n x logn lattice, the measurements can be simulated with poly(n) resources. 

Concerning the standard circuit model for quantum computation on pure states the points 1.-3. 
lead to: 

Theorem 19 (Simulating quantum circuits) If at every stage the N -partite state in a polynomial 
time circuit quantum computation has a MPS representation with D — 0{poly(N)), then the 
computation can be simulated classically with poly(N) resources. This is in particular true if 
every qubit-line in the circuit is crossed or involved by at most O(logN) two-qubit gates. 

For more details on the classical simulation of quantum circuits we refer to 1331 and 1341 for mea- 
surement based computation as well as for the relation between contracting tensor networks and the 
tree width of network graphs. The general problem of contracting tensor networks was shown to be 
jjP complete in 1351 . In 1361 a simple interpretation of measurement based quantum computation 
in terms of the valence bond formalism is provided, which led to new resources in 1371 . 
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A An open problem 

The main open problem we left open is the reasonability of assuming condition CI in the blocks 
of the canonical form. Since condition C2 is known to imply CI for sufficiently large Lo Q and 
we do have C2 in each block by construction, it only remains to show how big Lo can be. (We 
are neglecting complex eigenvalues of £ of modulus 1. The reason is that, by Thm[5] they can be 
avoided by grouping the spins into blocks or simply considering TV prime.) In the generic case this 
Lo can be taken Lo = [ 2 1 1 "j 3 j + 1. However there are cases in which Lo ~ 0(D 2 ). For instance 
Ao = 53 =1 I* + l)(*l> -^i = Our conjecture is that this is exactly the worst case. 

Conjecture 1 There exists a function f(D) such that for any sequence of D X D matrices Ai 
verifying CI for some Lo, Lo can be taken f(D). 

Conjecture 2 f(D) ~ 0(D 2 ). 

We have been able to prove both in a particular (but generic) case: 
Proposition 2 If Ao is invertible, then Lo can be taken D 2 . 
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Proof. Take an r and call S r C {0, . . . , d — l} r toa maximal set of indices for which 
{Ai 1 ■ ■ ■ -Ai r }(i 1 ,...,j r )£S is linearly independent. It is then sufficient to prove that the cardinalities 
#S r +i > #£V, whenever #S r < D 2 . If not, using that Aq is invertible, we can assume that 
Sr+i = S r x {0}. Now Ai y ■ ■ ■ A ir+2 = 

Ah C i2,...,'ir+2^'2 ' ' ' = 

(j2,---,Jr+l)6 S r 

E J2,...,ir+l fel,...,fcr 4, ...4, 42 

(j2,..-jr+i)ss r (ti,..,WeSr 

which implies that one can take 5V+2 = S r X {0} x {0} and hence #S r +2 = #5 r +i = #5 r . 
We can continue the reasoning and prove that, in fact, #Sk = #SV for every k > r; which is the 
desired contradiction. rj 

A particular case of this proposition occurs when one of the matrices in the canonical decom- 
position is hermitian. To see that one can group the spins in blocks of two. The new matrices are 
Aij — AiAj, and the canonical condition J]\ A\ — 1 reads £\ An = 1. Now, by doing a unitary 
operation in the physical system on can assume that Aq = y\ An = 1 (up to normalization) and 
therefore we are in the conditions of Proposition[2] 

The conjectures, if true, can be used to prove a couple of interesting results, one concerning the 
MPS representation of the W-state, and the other concerning the approximation by MPS of ground 
states of gapped hamiltonians 

A.l W state 

The Vl^-state can be written with matrices 2 x 2 in the OBC representation. However, this repre- 
sentation does not respect the symmetry of the state. What if we ask instead for representations of 
this kind? 

Let us start with a permutational invariant one, that is, one with diagonal matrices. By the 
permutational invariance of | Wn), it can be written as a sum X^^] 1 8n\x n )® N ■ Calling a n to a 
JV-fh root of S n and,4i = J2„ ct„(i\x n )\n)(n\, 

2 

\W N ) = J2 tr(A il ---A iN )\i 1 ---i N ). 

ill ■•■■'» = ! 

with the matrices Ai diagonal and of size (N + 1) x (N + 1). 

Of course we would be interested in a representation with smaller matrices. The surprising 
consequence of Conjecture [2] is that this is impossible even if we ask for a TI representation; that 
is, one with site-independent matrices (not necessarily diagonal): 

Corollary 1 If we assume Conjecture\2\and \Wn) = i N =i ^ r (^n ' ' ' ^»jv)I*i ' ' ' f or 

D X D matrices A„ then D>0 (iV 1/3 ). 

Proof. By Conjecture[2] we have condition CI by blocks with Lq = 0(D 2 ). Let us assume 
that y ^ 3(6 — l)(£o + 1), where (as usual) b is the number of different blocks in the canonical 
form. We can break the chain in two parts, each one having at least 3(fe — l)(Lo + 1) spins and 
decompose the state as \Wn) = X/f /3=i l^a^l^a^); where 

I$a,i3>= J2 ■ ■ ■ A iR \/3)\ii ■ ■ ■ i R ) 

i L ,...,i R 
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|*«,/3) = J2 (P\ A iR---Ai s \a)\i a +i---i N )- 

If we call Sj to the set of positions of the j-th block in the canonical form of \Wn) and S = 
USj, we claim that the sets {\^a,p)}( a ,i3)£S an d {\'&ct,p)}( a ,0)es are both linearly independent. 
This implies that the rank of the reduced density operator after tracing out the particles 1, . . . , R is 
— 53»=i Df, being Di x Di the size of the i-th block. 

We know that in the case of | Wn), this rank is 2, so we get exactly two blocks of size 1 X 1. It 
is trivial now to see that this is impossible. Therefore, < 3(6 — l)(£o + 1) and hence the result. 

So it only remains to prove the claim. For (a, j3) £ Sj let us take complex numbers c a ,p such 
that ~}2 a p c at f)\$ a ,i3) = ; which is exactly 



3=1 



E c <*ap){<*\ 

(c,/3)es 3 



By Lemma[4]the sum in j is direct, so each summand is 0. Finally, by condition CI by blocks 

J2( a ,p)eS- c a.p\P)i a \ = for every j and hence c a ,p = for every (a, (3) G Sj. \j 



A.2 Approximation of ground states by MPS 

As we saw in Sec(6] one of the big open questions in condensed matter theory is to mathematically 
explain the high accuracy of DMRG. Since DMRG can be seen as a variational method in the class 
of MPS, it is crucial to prove that any ground state of a gapped local Hamiltonian can be efficiently 
approximated by a MPS of low bond dimension D. Despite the important recent advances in this 
direction (see Sec(6^ the problem is still unsolved. An important step has been done by Hastigns, 
who has recently reduced this problem to the case in which the Hamiltonian is frustration free (5). 
This highlights the importance of the following dichotomy 

Theorem 20 (Dichotomy for the size of the MPS) IfConjecture&holdsandHN = Y^L X 

is a local TI Hamiltonian which is frustration free for every N, then the bond dimension D of any 

of its exact ground states, viewed as a MPS, is: 

( i) either independent of N 

(ii) or>0(Ni) 

Proof. Let us take the canonical decomposition of a ground state of H acting on TV 
particles 

d 

W)= tT(A il ---A iN )\ii---i N ). 

i-l,...,i N —l 

If (ii) does not hold, by Conjecture [2] TV > 3(6 — l)(£o + 1) + L (L the interaction length of 
H). Now, by Lemma[4] the products of the last TV — L matrices generate the whole space of block 
diagonal matrices. This immediately implies that the same matrices A, give us a ground state of 
H N i for any TV' > TV. That is, D is independent of TV. r- 1 
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